Fulde-Ferrell-Larkin-Ovchinnikov state in disordered 5-wave superconductors 



Qinghong Cui and Kun Yang 
National High Magnetic Field Laboratory and Department of Physics, 
Florida State University, Tallahassee, Florida 32306, USA 
(Dated: August 10, 2008) 

The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state is a superconducting state stabilized by a 
large Zeeman splitting between up- and down-spin electrons in a singlet superconductor. In the 
absence of disorder, the superconducting order parameter has a periodic spatial structure, with 
periodicity determined by the Zeeman splitting. Using the Bogoliubov-de Gennes (BdG) approach, 
we investigate the spatial profiles of the order parameters of FFLO states in a two-dimensional 
s-wave superconductors with nonmagnetic impurities. The FFLO state is found to survive under 
moderate disorder strength, and the order parameter structure remains approximately periodic. 
The actual structure of the order parameter depends on not only the Zeeman field, but also the 
disorder strength and in particular the specific disorder configuration. 



I. INTRODUCTION 

In the early 1960's, Fulde and Ferrellp and inde- 
pendently Larkin and Ovchinnikovp' proposed the 
possibility that superconducting states with periodic 
spatial variation of the superconducting order parameter 
would become stable, when a singlet superconductor 
is subject to a large Zeeman splitting. The Zeeman 
splitting could be due to either a strong external 
magnetic field, or an internal exchange field. Nowadays 
such states are collectively known as the Fulde- Ferrell- 
Larkin-Ovchinnikov (FFLO) state. Recently the FFLO 
state has received renewed interest due to experimental 
evidenc e of its existence in various superconduc- 
tors ^ ^ 5 I6I7I8I9I1QI11I12I13I14I15I16I17I1 8 19 20 21 22 23 24 25 26 27 

and the (thus far theoretical) possibility of its realization 
in trapped cold atom systems nuclear matter and 
in the core of neutron stars.^^ 

Soon after the original work of FFLO, it was found 
that the FFLO state is very sensitive to the presence 
of impurities, and eventually suppressed when t he disor- 
der strength reaches certain critical value .'^I'^^'^^^ Phys- 
ically, this is because the Cooper pairs in the FFLO state 
carry finite momenta (and the periodicity of the spatial 
structure is tied to the corresponding wavevectors); im- 
purities cause scattering between states with different 
momenta and thus tend to destroy the corresponding 
spatial structure. In these earlier studies however, the 
effects of impurities are taken into account by disorder- 
averaging various relevant physical quantities (like pair- 
ing susceptibility); after such averaging translation sym- 
metry is restored, and the superconducting order param- 
eter is again assumed to have a periodic structure. On 
the other hand in a specific realization of disorder, the 
system will adjust itself according to the impurity con- 
figuration,^^ and its order parameter structure will no 
longer be exactly periodic. We can thus expect (possi- 
bly quite large) spatial fluctuations of the order param- 
eter, which are neglected in these treatments. In fact 
such spatial fluctuations have been observed previously 
in strongly disordered superconductors in the absence of 
Zeeman splitting .'^^*^^ ' ^^ ' ^^ ' ^^ ' Thus a more complete un- 



derstanding of the FFLO state requires a more careful ex- 
amination of the disorder induced order parameter fluc- 
tuations. A related conceptual issue is how to distinguish 
the FFLO state from other competing states (including 
the Bardeen-Cooper-Schrieffer, or BCS state) in the pres- 
ence of impurities, where momentum is no longer a good 
quantum number. 

In this paper, we present a microscopic numerical 
study of a disordered, two-dimensional (2D) 5-wave su- 
perconductor subject to a Zeeman field. We use a nega- 
tive U Hubbard model and treat it within the Hartree- 
Fock (HF) and Bogoliubov-de Gennes (BdG) framework. 
We solve for the order parameter configuration specific 
to each disorder configuration, and develop empirical cri- 
teria to distinguish among the BCS, FFLO and normal 
phases. A phase diagram based on such criteria is ob- 
tained. 

This paper is organized as follows. In section II, we 
introduce the model Hamiltonian and present the mean- 
field treatment of our model. The numerical results of 
solutions of the order parameters are presented and ana- 
lyzed in section III. Some concluding remarks are offered 
in section IV. Throughout this paper we only consider 
zero temperature, although generalization to finite tem- 
perature is straightforward. 

II. MODEL AND MEAN-FIELD TREATMENT 

We begin with the disordered 5-wave superconductor 
on a 2D square lattice of finite size Nx x Ny^ described 
by the negative U Hubbard model: 

n = -t cj'^Cj^ + ^(^i + cr/l-/i)4^Ci^ 

<i,j>cr icr 

-Uj2ih,,-^-)inn-l), (1) 

where t is the nearest-neighbor hopping that will be set 
t = 1 from now on; a = ±1 is the spin index; Wi is 
the on-site random potential which is independently dis- 
tributed from —W/2 to W/2 uniformly; h is the Zeeman 
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field; U > is the attractive interaction; = Cj^Ci^r 
is the number operator, and i = {xi^yi) is the site posi- 
tion. Since the system size is finite, periodic boundary 
condition is imposed on the lattice. It is interesting to 
notice that the Hamiltonian ([T]) is symmetric under the 
particle-hole transformation 



(-ly 



(2) 



in combination with the transformations: fi — /i, h 
— h and Wi —Wi. This symmetry is respected by the 
mean-field approximation discussed below, and it is a 
very useful property that we use to test the accuracy of 
our numerical solutions. 

Within the HF-BdG approximation, the Hamiltonian 
reduces to a quadratic form 



•ArciiCii), 
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where i^ij = + (it^i -/i- [/^ni^)^ij; = (ct^Ci^) 

is the particle density and Suio- = nicr — 1/2; a = — <j; 
Ai = — [/(ci|Ci|) is the order parameter. To diagonalize 
this HF-BdG (or mean-field) Hamiltonian ([3|, we employ 
the Bogoliubov transformation 



(4) 



where and 7J are the quasiparticle operators. The 
amplitudes of the quasiparticles (u^cr^'^ia) satisfy the 
Bogoliubov-de Gennes (BdG) equations 



j 



] =eJ 



(5) 

where the eigenvalues Ej, >{). The self-consistency con- 
ditions are consequently expressed as 

Ai = -UY,[u\^v\^f(E,) - u\^vV;f{-E,% (6a) 

V 

ni. = + (6b) 

where j{E) is the Fermi distribution function. Self- 
consistent solutions of Ai and n\(j are achieved by it- 
eration of Eqs. (|6| from a randomly initialized configu- 
ration. In case different solutions are obtained from dif- 
ferent initial configurations, we choose the energetically 
favored one by comparing their corresponding energies. 
The energy of the system in the superconducting state is 
evaluated to be 



V icr i 



(7) 



At zero temperature the Fermi function f{E) reduces to 
a step function. 



III. NUMERICAL RESULTS 

In our numerical calculations, we take the chemical po- 
tential /i = — 0.8 to avoid the perfect nesting of the Fermi 
surface at half-filling (/i = 0), and a relatively strong pair- 
ing interaction strength U = 2.5 so that the FFLO state 
can be realized in a relative small system size. We start 
by considering the disorder-free limit of W = 0. In this 
case we can take advantage of the translation symmetry, 
and solve the problem in the thermodynamic limit. We 
find BCS order parameter A = 0.365 for = 0, and the 
FFLO state is stable for hc\ < h < hc2 with hd ~ 0.27 
and hc2 ^ 0.66. This field range is bigger than that of 
a 2D weak-coupling s-wave superconductor with circu- 
lar Fermi surface. We believe the discrepancy is due 
to the fact that our pairing strength is not weak, and 
lattice effects that result in a non-circular Fermi surface 
with some nesting that favors the FFLO state. The sys- 
tem turns normal (A = 0) for h > hc2- Throughout the 
range hd < h < hc2^ the FFLO order parameter takes 
a one-dimensional (ID) structure,!^ with a period that 
varies continuously with h. 

In the presence of W ^ we will study systems with finite 
size of rectangular shape, which is fixed to be = 32 
and Ny = 16 for the rest of the paper. Periodic boundary 
condition is imposed in both directions. Here, we make 
the lattice geometry to be rectangular instead of square 
because this allows us to extend the size of one of the 
directions (x-direction) at the expense of the other, and 
the ID structure of the FFLO state is not very sensi- 
tive to the size of the shorter direction. Having a finite 
systems size restricts the periodicity of the FFLO order 
parameter. At a given Zeeman field, if the lattice size is 
not commensurate with the corresponding period, solu- 
tions of the FFLO states will have to be modified to be 
commensurate with systems size, which results in some 
energy cost. 

In the absence of disorder {W = 0), FFLO state is 
obtained at this system size in the range from hd =0.27 
to hc2 = 0.55; the (slightly) reduced range refiects the 
finite size effect discussed above. The order parameters 
Ai at several Zeeman fields are plotted in the left column 
of Fig. [1] while the right column are the maps of the 
corresponding order parameters in Fourier space Ak, 

A.~EAie---. (8) 




Here, we see that the order parameter is of nearly co- 
sine form with period of 32 along x-axis at h = 0.27 
Jb)]. When the Zeeman field reaches h — 0.35 
Tic)], the period reduces to 16 (so that two peri- 
ods are accommodated), and it remains the same un- 
til the system becomes normal at = 0.56 [Fig. [ijd)]. 
In the cases of the FFLO state, we find two dominant 
Fourier components representing the corresponding sin- 
odal structure of the order parameter. 

To obtain the physical picture of the s-wave super- 
conductors subject to both disorder and Zeeman field. 
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FIG. 1: (color online) Order parameter at various h in the absence of disorder {W = 0) in both real space (left column) and 
momentum space (right column), (a) h = 0.0, BCS state; (b) h = 0.27, FFLO state with one period accommodated within 
the finite system size; (c) h = 0.35, FFLO state with two periods accommodated within the finite system size; (d) h = 0.56, 
normal state. Here, U = 2.5, /x = -0.8, A^^ = 32, Ny = 16. 
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we calculate the order parameters A with h = 0.24-0.70 
and W = 0.0-4.0. The disorder configurations are real- 
ized over 10 different samples. For comparison, we also 
calculate the order parameter Aq of each disorder con- 
figurations in the absence of Zeeman field, where Ander- 
son's theorem guarantees a BCS state within mean-field 
theory. In Figs. [2] and [3j we give some examples of the or- 
der parameter configurations with various Zeeman fields 
under weak and strong disorder strengths (W = 1.2 and 
4.0, respectively). Scrutinizing the spatial profiles of the 
order parameters individually leads to the following ob- 
servations. 

First of all, we find considerable sample specific spa- 
tial fluctuation in the order parameter in the BCS state 
[when /i = 0; see Figs. |2|a) andjsja)]. This is consistent 
with previous studies.^ On the other hand the order 
parameter A has a single (positive) sign despite the fluc- 
tuations, as a result of which the Fourier transform of A 
peaks sharply at k = 0; this is a hallmark of the BCS 
state. When h is turned on, there may be some sign 
changes in A; but as long as Ak=o remains the dominant 
Fourier component, we can still identify the state as BCS, 
and a typical example is Fig. (sjb). Due to the presence 
of a gap, the BCS state is unaffected by the Zeeman fleld 
for a flnite range of h. To determine the phase bound- 
ary across which the BCS state becomes unstable and 
yields to either the FFLO or normal state, we calculate 
the overlap between the order parameters in the presence 
and absence of Zeeman fleld, 

A = KAo|A)|/||Ao||^ (9) 

where ||A|| = ^/{A\A) is the norm of the order param- 
eter. The reason we focus on this quantity is that the 
overlap between the order parameters of BCS state and 
either FFLO state or normal state is exactly zero for clean 
system, because the FFLO order parameter is oscillatory 
with zero mean, and the normal state has A = 0. We 
expect A to decrease rapidly across the phase bound- 
ary where the BCS state destabilizes, and becomes al- 
most zero in the FFLO and normal states with disor- 
der. Numerical calculations of the sample averaged A 
are displayed in Fig. [4j There we can see a fairly clear 
phase boundary separating the BCS and non-BCS states. 
The boundary is sharper in the weak disorder regime, 
while becomes somewhat fuzzy with increasing disorder 
strength due to larger sample-to-sample fluctuations. 

We now turn to the more interesting FFLO state in 
the presence of disorder, examples of which are shown in 
Figs. [2|b), |2|c) and[2|d). Just like the disorder- free case 
{W = 0), the order parameters are oscillatory in these 
cases. But unlike the cases without disorder, the order 
parameters are no longer exactly periodic, and exhibit 
considerable sample speciflc spatial fluctuations similar 
to the BCS case with disorder. However the character- 
istic order parameter momenta are clearly visible in the 
Fourier space, as manifested by well deflned peaks in Ak 
away from the origin k = 0. Unlike the disorder- free 
cases however, these peaks have flnite width, reflecting 



the absence of exact periodicity. We also flnd that dis- 
order appears to influence the structure of the order pa- 
rameter, such as in Figs. [2|c) and[2|d), where the order 
parameters form 2D oscillatory structures while, with the 
same Zeeman fleld, only ID structures appear in the ab- 
sence of disorder. This is visible in real space but becomes 
particularly obvious in Fourier space, where Ak has four 
instead of two peaks. We thus flnd that the order pa- 
rameter structure depends not only on the Zeeman fleld, 
but also disorder strength and conflguration. 

Such disordered FFLO states are found at relatively 
weak disorder, up to W ^ 2.0. For stronger disorder a 
different pattern emerges as h increases. In this case the 
BCS state yields to a state in which the order parameter 
is zero or nearly zero almost everywhere, but there can 
be isolated regions (or islands) within which A is non- 
zero and quite random. In Fourier space Ak also shows 
a random form with no obvious structure, notably with 
Ak=o ^ 0. Examples of such states are Figs. |3jc) and 
[sjd). This pattern is similar to what happens when in- 
creasing the disorder to very strong strength at = 0, 
where previous work found that the order parameter is 
destroyed everywhere except in a few superconducting 
islands. The only difference here is that the order 
parameter tends to oscillate in sign in the presence of 
a Zeeman fleld. In Fourier space the order parameter 
shows random diffuse pattern with no clear structure in 
such states. Here we identify such states as the normal 
state, as in previous work .'^^*^^* This is because while the 
order parameter is non-zero in isolated islands, quantum 
fluctuation will easily destroy the phase coherence be- 
tween different islands, and as a result there is no phase 
coherence and stiffness in the thermodynamic limit. 

It is clear from the discussions above that these three 
competing states (BCS, FFLO and normal) are most eas- 
ily distinguished from the order parameter structures in 
Fourier space. To quantify their differences we introduce 
the following quantity: 

max /IIAkll, (10) 

where 77 = +1 if the maximum value of |Ak| (denoted by 
|Ak|max) is at k = and r] = —1 otherwise, and ||Ak|| = 
Y^^j^ l^kp- In the absence of disorder, B = 1 for the 
BCS state and B = -1/^2 for the FFLO state with ID 
cosine structure (with only two Fourier components; the 
magnitude of B reduces somewhat when more Fourier 
components are present). For the normal state, B = 0. 
In the presence of disorder, Ak is dominated by the k = 
component in the BCS state, and by several components 
with k 7^ in the FFLO state; we thus expect B to be 
positive for the BCS state, and negative for the FFLO 
state. For the normal state on the other hand, there 
will be no dominant component in Ak; as a result B is 
nearly zero for that case. To illustrate this identiflcation 
we have listed the values of B for the individual cases in 
Figs.[2]and|3l 

The disorder averaged B value is plotted in Fig. [5j 
Using the identiflcation above, we can clearly see these 
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FIG. 2: (color online) The same as Fig. [TJexcept W = 1.2 and (a) h = 0.0; (b) h = 0.32, A = 0.299, B = -0.601; (c) h = 0.36, 
A = 0.019, B = -0.447; (d) h = 0.40, A = 0.012, B = -0.548. See the text for the definitions of A and B. 
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FIG. 3: (color online) The same as Fig. [TJexcept W = 4.0 and (a) h = 0.0; (b) h = 0.34, A = 0.411, B = 0.511; (c) h = 0.42, 
A = 0.027, B = -0.195; (d) h = 0.64, A = 0.001, B = -0.099. See the text for the definitions of A and B. 
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h 

FIG. 4: (color online) The map of the overlap between A(/i) 
and Ao (/i = 0) , A = | ( Aq | A) | / | | Aq | | ^ , as a function of Zeeman 
field h and disorder strength W. We average over 10 samples 
of different disorder realizations for a given W. For the pure 
BCS state, the order parameter is unaffected by the Zeeman 
field and gives A = 1. For non-BCS state, A will decrease to 
zero rapidly. A rough boundary between the BCS state and 
non-BCS state is depicted by the (red) dashed line. 




0.0 0.27 0.65 h 



FIG. 6: Sketch of the phase diagram of disordered s-wave 
superconductors under Zeeman field, based on data presented 
in Figs. [4] and [5] The two critical Zeeman fields at = 
are based on our numerical results of the same U and /i as in 
Figs. |4] and [5] except that the lattice size is allowed to vary. 



three phases, and identify three phase boundaries: BCS 
to FFLO, BCS to Normal, and FFLO to normal. The 
first two boundaries agree with the BCS to non-BCS 
boundary obtained from Fig. |4] By combining data in 
these two figures as well as our knowledge of the posi- 
tions of hci and hc2 in the absence of disorder, we obtain 
a phase diagram of the system as presented in Fig. [6j 




FIG. 5: (color online) The map of B = '/7I Ak|max/| | Ak| | as 
a function of Zeeman field h and disorder strength VK, where 
Ak is the Fourier transform of Ai, 77 takes the value of 1 when 
the maximum value | Ak|max is at k = 0, and of —1 otherwise. 
The disorder realizations are the same as in Fig. [4] In the 
absence of disorder, this quantity equals to 1 for the BCS 
state, —l/y/2 for the FFLO state of ID structure while —1/2 
for 2D structure with simple cosine functions for the order 
parameter, and zero for the normal state. The (red) dashed 
lines sketch the (rough) phase boundaries among these three 
phases with disorder. 



IV. CONCLUDING REMARKS 

In this paper we have presented a microscopic study 
(based on mean-field theory) of a 2D s-wave supercon- 
ductor subject to both disorder potential and a Zeeman 
field, and showed that the FFLO state can survive mod- 
erate disorder strength. Our most important results are 
presented in the form of a zero temperature disorder- 
Zeeman field phase diagram (Fig. [6|, as well as the de- 
tailed analysis of the superconducting order parameters 
in different phases that leads to the determination of this 
phase diagram (Figs. [4] and [5|. We should note that 
within mean- field theory we always obtain a BCS super- 
conducting solution for zero Zeeman field (/i = 0); how- 
ever quantum fluctuations that are left out in our study 
gets enhanced with increasing disorder strength and 
will eventually drive the system to an Anderson insulat- 
ing phase. Thus the BCS-Normal phase boundary will 
eventually bend left as W increases, and intersect the 
h = axis. Since we do not expect metallic behavior 
in 2D, the entire normal phase is expected to be an An- 
derson insulator. We also note that the superconducting 
islands we found in the "normal" phase (which is really 
insulating) have also been observed in numerical stud- 
ies of superconductor-insulator transition in disordered 
superconductors.'^ 

Our results are in qualitative agreement with previous 
studies^^ which suggested that the the FFLO state can 
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survive as long as the electron scattering rate is compa- 
rable or smaller than the BCS gap. In our study we find 
the critical disorder strength Wc ^ 2.0; this corresponds 
to a scattering rate (within effective mass approxima- 
tion) I/tc ^ W^/{2m) ^ 0.167 (we set h = 1), which 
is indeed comparable the the BCS gap A = 0.365. More 
importantly however, our study reveals that the FFLO 
state is characterized by approximate^ instead of exact 
periodic variation of the superconducting order parame- 
ter in real space with zero average; this is manifested in 
Fourier space by peaks in the order parameter at non- 
zero wavevectors with finite width. 

In an earlier work,'^one of us (along with D. F. Agter- 
berg) proposed to use the Josephson effect between a 
BCS superconductor and an FFLO superconductor to 
identify the latter and measure its order parameter struc- 
ture. It was shown thalP^ there is no Josephson effect in 
this case when there is no magnetic fiux going through 
the junction, because there is no k = component of 
the order parameter in the FFLO state. This conclusion 
remains to be true in the presence of disorder. Ref.(4Tfur- 
ther predicts that the Josephson effect can be recovered 
when there is an appropriate amount of magnetic fiux 
going through the junction, which gives rise to a momen- 
tum (or wavevector) to the Josephson hopping matrix 
element that matches momentum of one of the Fourier 
components of the FFLO order parameter; this allows for 
a direct measurement of the order parameter momenta in 
the absence of disorder. In the presence of disorder this 
method should still work as long as the disorder strength 
is sufiiciently weak and the size of the junction is rela- 
tively small. However for stronger disorder and/or very 



large junction this will no longer work because the weight 
carried by a single Fourier component of the order param- 
eter will be negligible in the disordered FFLO state. 

It would be highly desirable to go beyond mean-field 
theory and take into account effects of quantum fluctua- 
tions of the order parameter. This is perhaps most eas- 
ily done in one- or quasi-one-dimensional systems, where 
the powerful method based on bosonization has already 
yielded a number of exact results in the absence of dis- 
order. Disorder can be treated using the formalism ad- 
vanced by Giamarchy and Schulz.^^ We note that in 2- 
or 3-dimensional systems, initial attempts to go beyond 
mean-field theory have been made by treating fiuctua- 
tions within RPA^ and fiuctuation exchang^^ approxi- 
mations. 

Among existing candidates of FFLO superconductors, 
CeCoIns appears to be the most promising. This is 
a (i-wave superconductor. The d-wave FFLO state is dif- 
feren t from g -wave FFLO state in a number of important 
^g^yg"34|49i50j would thus be very interesting to use the 
methods of the current paper to study (i-wave supercon- 
dutors in the presence of both disorder and Zeeman field, 
and compare with existing theory. 
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